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The decay and annihilation cross-sections ol dark matter particles may depend on the value of 
a chameleonic scalar field that both evolves cosmologically and takes different values depending on 
the local matter density. This possibility introduces a separation between the physics relevant for 
freeze-out and that responsible for dynamics and detection in the late universe. We investigate how 
such dark sector interactions might be implemented in a particle physics Lagrangian and consider 
how current and upcoming observations and experiments bound such dark matter candidates. A 
specific model allows for an increase in the annihilation cross-section by a factor of 10 between 
freeze-out and today, while different choices of parameters allow for scattering cross-sections near 
the astrophysical bounds. 



I. INTRODUCTION 



The particle physics properties of dark matter arc important for three distinct aspects of its behavior: they determine 
how the initial abundance of dark matter arose, they govern how the dark matter distribution evolves and influences 
structure formation, and they delineate the possible ways in which dark matter may be detected. Of course, these 
three roles are not typically independent, since they all depend on the prescribed interactions between the dark 
matter particles themselves and also between dark matter and the Standard Model. These connections often provide 
a powerful motivation for particular dark matter candidates - for example, the freeze-out abundance of WIMPs points 
to new physics at the weak scale, which in turn leads to an attractive connection between dark matter and proposed 
solutions to the hierarchy problem, such as weak-scale supersymmetry. 

The idea that dark matter could have interactions of astrophysically interesting magnitude has received a good amount 
of attention [IHZ], motivated in part by purported discrepancies between the standard ACDM model and observations 
of structure on small scales (as described in [5], for example). One obstacle to simple implementations of this idea is 
that the required cross-section for a thermal relic to obtain the right relic abundance is close to the weak scale, far 
too small to be relevant to dynamics in the late universe. In this paper we explore the idea that the dark matter 
cross-section might be much larger now than it was at freeze-out, due to the evolution of a background field. 

In a cosmological context, the evolution of background fields can assert a significant influence on the properties of 
dark matter as a function of spatial location or cosmic epoch [9HT8]. A straightforward way to achieve such effects 
is to invoke a light scalar field that interacts with dark matter and/or ordinary matter as well as through its own 
potential, and whose expectation value feeds into the dark-matter properties. A popular scenario along these lines is 
the "chameleon mechanism," which acts to screen light, cosmologically relevant degrees of freedom to protect them 
from precision local tests of gravity [I9ll23j . 

In this paper we investigate dark matter that interacts through a gauge symmetry with a coupling constant that 
depends on a chameleon-like scalar field. (The effects of chameleon vector bosons on laboratory experiments were 
considered in 24 .) Just as the properties of a cosmologically relevant scalar can be drastically modified in the 
presence of local density inhomogeneities or after evolving over cosmic time, so the interactions of dark matter may 
be modified. We are able to find a model in which the late-time interaction strength is considerably higher than that 
at freeze-out - although admittedly, this behavior does not seem generic (without violating experimental constraints 
of one form or another) . 
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We begin by re-examining the conventional story of dark matter freeze-out according to the Boltzmann equation, 
but with the additional ingredient that the dark matter properties are evolving with time. We then look at specific 
models featuring a Dirac dark matter particle and a U(l) gauge symmetry that is spontaneously broken, along with 
a chameleon scalar field. We study the cosmological evolution of this coupled system, and calculate the dark matter 
properties, including annihilation and scattering cross-sections. Finally we exhibit numerical solutions to a specific 
model, showing that the annihilation cross-section can increase substantially during cosmic evolution. 



II. THE GENERAL PICTURE: EVOLVING DARK MATTER IN THE EARLY UNIVERSE 



Before discussing specific models, let us first consider how the usual story of dark matter freeze-out might be modified 
if the annihilation cross-section depends on the dynamics of another field. In the next section, we will explore La- 
grangians that couple the dark matter to a scalar field that affects its interaction cross-sections. For simplicity we work 
in a fiat Friedmann, Robertson- Walker (FRW) universe, described by the metric ds 2 = —dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ), 
with scale factor a(t). 

The decoupling of dark matter takes place in the early universe in the radiation-dominated regime, in which particles 
with masses m <C T are the dominant component of the cosmic energy budget. To a good approximation, we may 
therefore ignore contributions from non-relativistic species in thermal equilibrium with the radiation and approximate 
the energy density as 



7T 
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and the entropy density as 



where, as usual, 
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and gi is the number of internal degrees of freedom for particle species i. 

For T > 300 GeV, g* s = g* = 106.75, which includes all particles in the Standard Model. When 100 MeV > T > 
1 MeV, the electron and positron are relativistic and so g*s — g* — 10.75. At the temperature of the CMB today, 
T = 2.725 K, g* St0 — 3.91, and g* fi — 3.36. 

Consider a dark sector that was in thermal equilibrium with the visible sector at some very high temperature scale, 
below which they decouple effectively enough to consider each sector separately to be in equilibrium. The visible 
sector is at temperature T with entropy density s(T), while the dark sector is at temperature Td with entropy density 
Sd{Td)- The expansion of the universe is governed by both sectors with 

9T{T)=g^T)+gt{T d ){^ , (5) 

but quantities in the dark sector (for example, the dark matter annihilation cross section and number density) are 
determined by Td [25] , 

Since the entropy in each sector is conserved independently, the assumption that the two sectors were in equilibrium 
at some unification scale at time t u allows us to express the dark bath temperature in terms of the visible bath 
temperature at some later time t via 
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All Standard Model particles contribute at t u to give g*s(tu) — 106.75, and all dark particles contribute to gfg(t u ). 
In what follows, we will use the temperature of the visible sector and convert to T as needed. For convenience we 
write 

m = TM = ( g*s(t) gf s (t u ) \ 1/3 ^ (7) 



T(t) W s (t) 9*s(t u ) 



The success of Big Bang Nucleosynthesis (BBN) and the structure of the Cosmic Microwave Background (CMB) 
power spectrum place tight bounds on any new relativistic degrees of freedom in the dark sector. The limit on the 
effective number of light neutrino species is N v — 3.24 ± 1.2 at the 95% confidence level [26], which gives 

S*£ 4 0bbn) = I x 2 x (N v - 3) < 2.52 (95% confidence) (8) 
8 

for 3 light SM neutrino species [37]. The 5-year WMAP data [25] also bounds the number of neutrino species by 
N u = 4.4 ± 1.5 at the 65% confidence level, and the 7-year WMAP data [5S] places a tighter lower limit of N v > 2.7 
at the 95% confidence level. 



A. The Boltzmann Equation 



Let us assume the dark matter ip is a stable particle that annihilates with a thermalized annihilation cross section 
(av) . The general Boltzmann equation governing the number density n of a particle of mass m is 

h+3Hn+ (av) (n 2 - n| Q ) = , (9) 

where H is the Hubble parameter 



H = a = \jr GpR = i^^ T (10) 

and riEQ is the equilibrium number density 

where K 2 is the modified Bessel function of the second kind of order two. Generalizing the traditional treatment, we 
allow for the possibility that the mass of the dark matter m^,((/)) is a function of a real scalar chameleon field <p and 
denote </>-dependent masses and couplings with a tilde. 

It is convenient to scale out the effects of the expansion of the universe by defining 

Y = ^ (12) 
s 

{n^,{x) and Y(x) are taken to be independent of 4>) and to use a new independent variable, related to the cosmic time 
t through 

x{t) = ^ , (13) 

where mj is some constant mass scale. In the usual derivation, mj is chosen to coincide with the dark matter mass; 
however, since our dark matter has varying mass, we use this constant parameter instead. Defining 



47r 3 G m T 
allows us to write 

dx rriT 



(14) 
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which can be used to rewrite the Boltzmann equation for the dark matter as 

B 



Y>(x) + Yi Q )=0. (16) 



x 



Here a prime denotes a derivative with respect to x, and 



5=M^^<, (17) 



45 y^tot 

which may depend implicitly on <j) in our model via a ^-dependence in the cross section. Note that, in terms of these 
new variables, the equilibrium term is 

45g / ( x m^((f>) 



with g = 2 for Dirac dark matter. 

It remains, at this level, to specify Y(xi), the initial condition for Y. We consider A = Y — Yeq, the departure from 
equilibrium |30j . which obeys 

A' = -Y^ Q - ^A(2F EQ + A) . (19) 

At early times (1 < x <^ Xf), Y tracks Yeq extremely closely such that A and |A'| are small. Note that in the 
non-relativistic approximation, T <C m^, (</>), 

Yeq ~ x^^e-^/Oi^W/^T) > (20) 
and so ^eq/^eq ~ — w^t^/^T Emd A' ss 0. Thus, the required initial condition is 

r <-» = ^<-» + 1g^' < 21 > 

where B((f>i) and m^,(<^>j) are evaluated at the initial value <f>i = (j){ x i)- 

After the freeze-out value Xf, Y(x) will asymptotically approach a constant value Y^. The energy density of non- 
rclativistic dark matter today is then 

Pa = m^(4>o)n^(x ) = m^(0 o )^ooSo 
2?r 2 

= m lp ((j) )Y oo —g^ Sfi T^. (22) 

Having generalized the usual treatment of dark matter as a fluid to the case in which there is a chameleon field 
determining the dark matter properties, we now turn to specific examples of particle physics models in which these 
phenomena might arise. 



III. GAUGED DARK MATTER 



Consider dark matter to consist of a Dirac fcrmion tp, charged under a dark U(l) gauge group with gauge boson 
A^, and a dark Higgs field $ that spontaneously breaks the U(l). We also introduce a chameleon field </> that is a 
real scalar field with properties that depend on the dark matter energy density. The chameleon couples to the other 
particles in the dark sector by entering into the dark matter mass m^, ((/>), the U(l) coupling /(</>), and other couplings 
described below. 
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A. A Toy Model for Varying Coupling 



As a first step, let us consider the QED Lagrangian with a real scalar held (f>, but in which we allow the coupling 
constant e to vary as a function of spacetime |31j . Specifically, it can vary as a function of <p. Let us write the new 
coupling as /(</>). Thus, 

£qed0 = -\d^d^ - V{4>) - -^—F^F^ + i$H - - ^ 7 M V^ , (23) 

1 VP) 

where — d^Ay — dyA^. Making the redefinition A^ — \ f{4>)A^ : we obtain 

£qed = -\d^(t> - V(cf>) + ii>H - - fWftfliA,, - -I- \dP(fA v ) - d'iUn} 2 • (24) 

z 4/ 

Both Lagrangians are equivalent, but now the gauge transformation reads 

f(<p)A^ -)■ f{cp)A^ + d^cu (25) 
ip -> e^i/) (26) 
ip -> e +iu ip . (27) 

If we can neglect factors of (<9 M / / /) compared to all other mass scales in the theory (except the Planck mass), then 
the Lagrangian simplifies to the approximately gauge-invariant form 

£ QED « ~\d^4> - V(<t>) - \f^F^ + ~ - f{^M (28) 

with U(l) current 

f{x) = f(cp)^ . (29) 



B. The Cosmological Equations of Motion 



We now include gravity and a complex dark Higgs field $ to break the U(l) symmetry and give the dark gauge field 
a mass. We allow for a varying dark matter mass by using the effective mass parameter rh^p(4>), and in the spirit of 
effective field theory, we also allow all couplings (not just the U(l) coupling /(</>)) to depend on (f>. 



Neglecting factors of {d^f/f), the action is then 



n 



16ttG 
1 

~ 4 



V{<j>) ~ - V ($) 

F^F^ + ~ m^cp)^ - A^(0)($ + $t)^ 



(30) 



where the gauge covariant derivative is = V M + if ((f)) A^. The equations of motion for the fields then follow as 

^0-m^(0)-A^(0)($ + $ t ))^ = O (31) 

- v'(4>) + ffyttyw - f{^H - a;(£)($ + * f )v^ = o , (32) 

where a prime denotes differentiation with respect to <fi. Let us assume that the universe is dark-charge symmetric, so 
the average charge current density is negligible compared to the dark matter number density (see (37) below). Thus, 
the term proportional to /'// should be small compared to the one containing m'/m, given that /'// ~ m'/m to 
within a few orders of magnitude - a condition we will enforce later. We may write this last equation as 



□ 



- V'(cf>) + m'44>)^ - A;(0)($ + $t)^, ~ o 



(33) 



G 



We will arrange for the dark Higgs to have a sufficiently large mass that its perturbations are negligible and simply 
replace $ by ($) in the equations of motion. The VEV generates an additional mass term for ip, but we can redefine 
rh^{tp) to absorb this term. We then have 

(0-m^O)) ip&O (34a) 
Utp-V'{tp)+fh'Jtp)ipipnQ . (34b) 



We calculate the energy-momentum tensor for ip by varying the action with respect to the metric. Taking care to 
correctly handle the non-trivial metric-dependence of the covariant derivative |32j , we have 



» v 2 

where we have integrated by parts and used the field equation of motion. Taking the trace, we obtain 

i r_ . -<-r 



(35) 



1 fJLV 



Ipflp-lpflp - f(cp)lp£lp 



fh^(<p)ipip , 



(36) 



where, again, we have used the Dirac equation for ip and ip to obtain the last line. If we model the dark matter as 
non-relativistic dust, its pressure is zero and so the trace of the stress-tensor is approximately given by — p^. Thus, 



p$ = ~rh^(tp)ipip 



As a final step in this section, we use this result to rewrite the (^-equation of motion ( 34b ) as 

□0 - K' ff (0) = o , 

where the effective potential is 

Vcs = V(<f>) + m^{4>)n^ 

= V(cP) + rh^<p)Y(x)^g* s (^f ■ 



(37) 
(38) 



(39) 



IV. CHAMELEON BEHAVIOR 



With a complete model in place, we now turn to a detailed investigation of the dynamics. We first examine the 
chameleon field, which is central to the effect we seek. Assuming that <p is homogeneous and isotropic, so that we can 
neglect spatial derivatives in □</>, the equation of motion becomes 



4> + 3H<P + V'{4>) + = ■ 

It is convenient for seeking numerical solutions to work with a dimensionless variable 



(40) 



P = 



rriT 



and to use x as our independent variable. The equation of motion becomes 



P"(x) + -P'{x) + 



dV 



x ™yff' ot 



(ji—PrriT 



2ir 2 b 2 g*s drh^ 
45a; gl ot dtp 



Y(x) = 



(41) 



(42) 



(j)—PTYlT 



We choose the initial conditions for tp to begin at the minimum of its effective potential and to move with same 
initial velocity as the changing minimum. The minimum tp m i n solves the equation V^g{tp m i n ) = 0, so one of the initial 
conditions for this equation can be obtained by evaluating this expression at Xi, using the relevant value for Y(xi) 
from (21). Furthermore, since tp m - m is a function of x, the initial velocity is found simply by taking a derivative and 



using the Boltzmann equation to obtain the relevant value for Y'(xi). 
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A. Exponential Potentials 



Our goal here is to work out a single example model that exhibits the effects we are investigating, while at the same 
time remaining compatible with experimental constraints. For simplicity we will choose exponential functions, which 
also have the nice feature that observables approach a fixed asymptotic value at late times. 

With these comments in mind, we therefore choose the form of the effective potential and U(l) coupling to be 

V(4>) = A 4 e~^ mi (43a) 
m*W=m^ (l-A 2e -^™ 2 ) (43b) 

f(q>) = e(l + A 3 e^ m A~ 3 , (43c) 



where A and my, are constants with dimensions of mass, and e and A2, A3 > are dimensionless. The largest energies 
of the particles in our theory are of order my for non-relativistic dark matter, since all other particles should be lighter 
than the dark matter to allow for annihilation. We, therefore, require m,2, m 3 3> my, to suppress higher-dimensional 
operators involving derivatives of my and / when we expand the action. Additionally, we need nil > A to suppress 
higher-dimensional operators in the self-couplings of </>. 



The effective potential in ( 39 ) is now 



V eS (<j>) = A 4 e-^ + my,(l - A 2 e-^)Y(x) 2 ^g irS {^ff , (44) 
possessing a critical point at 



mim 2 /. m 1 m^m^Y'. 

ln i A 2- . (45) 



rmm I ^ a 4 

ni2 — mi \ m2 A* 

which is real and finite. In order to generate a mass for the excitations of 0, we require this critical point to be a 
minimum, which holds for 

m 2 > mi . (46) 

The minimum moves with a speed 

mima ^_3\ (47) 



dx m.2 — mi \ dx x 

which is positive (^> m in increases with x). Finally, we identify the initial conditions for </>: 

> mim 2 ( mim lp m^Y(x l )\ 



(49) 



TO2 — mi \ m 2 A 4 
(a;,) = - mim2 [ r(x i ) < ^-(a; i ) - — 



dx m2 — mi \ dx Xi 

mim 2 ( 3 my, Ieq^) + {xfm^)/(4B^m T ) 



m 2 -mi V^i £mr F E Q(xi) + (x^fh^)/(2B^m T ) 
In order to ensure my, > 0, we require 



> m 2 ln(A 2 ) (50) 



for all (/> relevant for our calculation. 



B. An Attractor Solution 



A particularly interesting and simple possible evolution for the chameleon field is for it to begin at the minimum 
of the effective potential and then to adiabatically track this minimum as it evolves cosmologically. This attractor 



solution [53] is achieved if the physical mass of the chameleon satisfies 



Ph = VKWmin) » H . (51) 



If (51 1 holds during radiation dominance, when 



H R = n YVW t ^~\ (52) 



then we can avoid solving the coupled differential equations (16) and (40) and simply use the expression for </> m i n for 



the evolution of (fr. Similarly, if (51 ) holds during matter domination, when 

3/2 



Hm — H { 



>(W ■ <»■> 



then we can easily determine (fro, the value of <fr today, which is needed to calculate the values of the ^-dependent 
parameters today. 

Under the approximation that TO2 3> mi, 

^nU^^^Y^ggx-*'*. (54) 
\ 45 m 1 m 2 / 

It follows that Hr decreases more rapidly than m^ p h with time, whereas during matter domination, Hm and m^ph 
have the same ^-dependence. We shall verify later that these attractor solutions exist by numerically solving all the 
relevant equations of motion. 



V. PARTICLE PHYSICS INTERACTIONS AND CONSTRAINTS 



In the adiabatic regime described above, we now have all the ingredients necessary to understand the cosmological 
evolutions of the fields. We next turn to the particle physics phenomenology of the model. To do this, we rewrite the 



action ( 30 ) without gravity to give the Lagrangian 

1 

-F^F^ + i$$*l> - rh^((fr)4,tfr - f(<fr)i>M - A^)($ + , (55) 



C « - td^tf, - V((fr) - (D^(D^) - Vo($) 



with D„ = d fl + if{</>)A„ 



A. Breaking the Dark (7(1) Symmetry 

The potential of the dark Higgs field $ is chosen so that this field acquires a vacuum expectation value (VEV) 

(0|$(x)|0) = . (56) 

Decomposing $ into two real scalar fields via 

*(a:) = 4=(v + h{x))e- lx ^' v , (57) 
we can then use unitary gauge x( x ) = to rewrite the kinetic term for $ as 

- (D^D^ = —d^hd^h - \f 2 {cfr){v + hfA^A^ . (58) 
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Thus, the Goldstone boson is eaten to give the dark U(l) gauge boson a mass M A {<f) = /(<^)v. The Yukawa term 
generates a contribution to the mass of ip, but since already has a Dirac mass, we need not rely on the dark Higgs 
to be the sole source of the ijj mass. We, therefore, absorb the dark Higgs contribution into the definition of rh,p and 
retain the freedom to choose this mass scale and the coupling X^(4>) separately. 



A typical choice for the pure dark Higgs potential Vo(<fi) is 

F ($) - \\ h (4>) 
which, when expanded about the VEV, yields 



~] 2 



$t$_ _ v 2 

2 



V (h) = - A \ h ^)v 2 h 2 + - A \ h {^h z + ^\ h (cj>)h 4 



(59) 



(60) 



The mass of the physical dark Higgs particle h is therefore 



and we see that the masses of the A and h fields are then related by 

2 



M A {4>) = f(4>\ 



rhh{<t>) 



(61) 



(62) 



Since the relative sizes of f(4>) and Xh(<P) are unrestricted, in principle the relative masses of A and h are not fixed. 
However, in order to simplify the analysis, we will impose the hierarchy rhh((f)) > 2M A {<j)) for all relevant <fi so that h 
has a tree-level decay channel to A. 



Our Lagrangian at this stage is then 

£ = - Idu 
2 M 

1 



1 



l r_ r. 



hA^A,, 



2/(0)M A (0) hA^A^ - V2X^(4>)h^ . 



(63) 



What remains is to incorporate the fact that (f> is adiabatically tracking the minimum of its effective potential. To 
achieve this, we expand 4>( x ) — 4>c(t) +vi x ) around its classical value, and recall that m 2 and 7713 are sufficiently large 
to suppress non-relevant terms of 0(rj) or higher. The Lagrangian (63) then becomes 



1 fr 



\ h (<t>c)+0{r,) 
M 2 A ^ C )+0( V ) 



v 2 h 2 - 
A^A, 



V{<t>c) 

v/1 3 



1 

F 2 
1 

24 



V"{cj> c )r? + 0{r?) 

|a a (0 c ) + 0{n) 



f{4> c )+0{rj) ^A^ 



[m^ c ) + fh^(cj> c )r} + 0(rj 2 )] ^ - ^2 [\^(<f> c ) + 0(rj) 

1 



- \2f\<j> c ) + 0{r,) 



h 2 A"A, 



2f(4> c )M A (cj> c ) + 0(r)) 



hA^A,, 



(64) 



B. The Dark Matter Annihilation Cross Section 



Our central goal is to understand how the dependence of dark matter cross-sections on the chameleon field changes 
the standard dark matter creation, evolution, and detection story. To this end, we next turn to the calculation of the 
dark matter annihilation cross-section. The relevant Feynman rules can be found in the appendix. 
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iMi 






iM 2 = 








FIG. 1: Tree-level ijj annihilation diagrams. The massive vector boson A is a wavy line, and the scalar h is a dashed 
line. Annihilations to A + A and h + h via ^-exchange and annihilations to final-state rj particles are suppressed by 

large-mass factors. 



We assume that the dark matter is the heaviest particle in the dark sector, such that rh^ 3> rhh, Ma- Then, the 
lowest order, tree-level processes for 2 — > 2 dark matter annihilation are shown in Fig. [TJ and their amplitudes are 



v / 2/(0 c )A^(0 c )Ma('/' c )A/ 1 (pi + p 2 )g^ ui 



(65) 



M 2 = iv 2 



2A$(0 C ) (A^(pi - k[) + A^( Pl - k' 2 )) 



3 - 

—T=X^{(j) c )X h {(l) c )\l/^h{Pl +P2) 

V2 



"i 



(66) 



M 3 = ie Vv v 2 [\/2/(<£ c )A^ c ) (A^(pi - k[)j v + l v ^(px - k[)) 

+ Pi^MAi^h^ipx +p 2 )] ui • (67) 
Working in the center-of-mass frame and in the non-relativistic limit we then obtain 

» 1fi ~2 lrk \ ( 68 ) 

167rm^(0 c ) 
15A^ C ) 2 

~ 7^; — I o ; , x « (69) 

128^m2(0 c ) v 7 

/ 2 (0 C )A^ C ) A 2 (0 C ) 

a-iV « = — = — — , (70) 

8ttM2(0 c ) 8™ 2 v ; 

where v is the relative velocity of the incoming particles. 
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(a) Annihilation (b) Scattering 

FIG. 2: Ladder diagrams for dark matter annihilation and scattering. 



The interaction of main interest is between the dark matter and the gauge boson mediator. The dark Higgs's primary 
role is to break the U(l) symmetry to give the mediator a mass, and most of its particle interactions can be neglected. 
The contribution &i is p-wave and thus sub-dominant to the other processes, which are s-wave. Also, the diagrams 
involving exchanges of h in Mi and M2 do not significantly contribute. Thus, the dark Higgs has the opportunity 
to influence dark matter annihilations only via 03. Let us insist that the Yukawa coupling Aw, is small enough (recall 
that the dark matter does not rely on this coupling to obtain a mass) such that 03 can be safely ignored. In this case 
only ci remains and, since it is an s-wave cross section, it is a simple task to carry out the thermal averaging required 
in the Boltzmann equation. Note, however, that if thermal averaging is needed (following Ref. [34]), we must use the 
dark sector temperature Td in the expression 



1 



9 



(n* Q (T d W 2 ( 2 -) 



ds \fsT c iK\ 



T,, 



(71) 



C. Corrections to the Cross Section 



We are interested in non-relativistic dark matter, for which the relative velocities are much less than the speed of 
light. It is well-known that for sufficiently low velocities, non-perturbative effects can have a large impact on the 
annihilation and scattering cross sections; and ladder diagrams, such as the ones shown in Fig. [2j must be included 
in the calculation. 



1. The Annihilation Cross Section 



In the case of annihilation, performing this summation is equivalent to solving the Schrodinger equation in quantum 
mechanical scattering theory [35] . This yields the so-called "Sommcrfcld enhancement" [36] of the annihilation cross- 
section (for detailed reviews in the context of dark matter, see, for example, [37-39 ). We consider the annihilation 
cross section <ro for a point-like interaction near r = in perturbative field theory. For small velocities, the attractive 
Yukawa potential 

V(r) = - a e -^ Ar , (72) 
r 

where a — f 2 (4> c )/Ai{, distorts the wave function at the origin and cannot be ignored. Including the potential will 
enhance the annihilation cross section to a — aoSk by the Sommerfeld enhancement factor Sk- Let us define the 
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dimensionless parameters 

e v = | (73) 

ex = ^ , (74) 

where v is the velocity of each annihilating particle in the center-of-mass frame. In the case of a massless gauge boson 
with a Coulomb potential, it is possible to solve the Schrodingcr equation analytically to obtain the Sommerfeld 
enhancement. 

For a massive gauge boson, the situation is more complicated, since the attractive potential has a finite range that 
limits the enhancement from being arbitrarily large for very low velocities. In the regime ca -C e 2 , we recover the 
Coulomb case. At the crossover point e v ~ ea (or equivalently rh^v ~ Ma), the de Broglie wavelength of the dark 
matter becomes comparable to the range of the interaction. At lower velocities with €a ~> e^, the Yukawa potential 
cannot be ignored. As v — > 0, the de Broglie wavelength increases to a value larger than the interaction range, and 
thus the enhancement saturates at 

S *~i~^. (75) 
ca m^, 

Furthermore, zero-energy bound states may form for certain values of ex, giving resonance regions with larger en- 
hancements ~ ^a/^v until they are cut off by finite width effects. In the early universe, freeze-out typically occurs 
at velocities Vf <~ 0.3, so that e v > 1 and the Sommerfeld enhancement can be ignored. Note that there are no 
enhancements for ca > 1- 

To find the thermally averaged cross section, taking into account the Sommerfeld enhancement, we integrate Sk using 
a Maxwellian distribution 

F(v) = ^—v 2 e-" 2 /* 2 , (76) 

where v is the characteristic velocity of the astrophysical system of interest. Thus, 

(av) = (<7u) s _ waV c (Sk) (77) 

/"CO 

(S k ) = / dv F(v)S k . (78) 
Jo 



For the purposes of this paper, we choose to work in the > 1 regime. This has two consequences. Practically, the 
calculation becomes much simpler, since we need not worry about the Sommerfeld enhancement at all. In addition, 
by deemphasizing the Sommerfeld enhancement, we clarify the extent to which the novel effects developed in this 
paper can alone increase the cross section over time in areas of parameter space that the Sommerfeld enhancement 
cannot reach. 



2. The Scattering Cross Section 



To find the scattering cross section, we can use non-relativistic quantum mechanics and sum over partial waves. The 
total cross section is 

A 00 

CT = ^E( 2Z + 1 ) sin2 ^ . ( 79 ) 

(=0 

although a more useful quantity to compare to observational constraints is the transfer cross section 

da 



(T tr = J dfl (1 - COSl,^ 

Air 



= T2 5Z[(2/ + l)sin 2 ^-2(; + l)sin^sin^ +1 cos( ( 5 i+1 -<*,)] , (80) 
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FIG. 3: Tree-level h decay. 



which controls the rate at which energy is transferred between colliding particles. Following [40j , analytic estimates 
for the cross section are 

4-7T 

<?=^^(l + L) 2 (81) 

^ = 4^(1 + L), (82) 

where L = /iu rc i& ma x is the largest angular momentum needed to describe the interaction between two particles of 
reduced mass \x = that travel with a relative velocity v Te i and maximum relevant impact parameter 6 max . Note 

that these estimates are only valid for L > 1. We estimate the impact parameter by solving 

= f^e-***~ . (83) 

"max 

If we work in the ca > 1 regime to avoid Sommerfeld enhancements, then we will also tend to avoid enhancements to 
the scattering cross section and can expect to be working in the Born limit. Simply taking the non-relativistic limit 
of the pcrturbative cross section gives 

/ 4 (0eK(0c) faA . 

a = ^— . (84) 

&TtM%{<t> c ) 

Assuming that dark matter self- interactions is not needed to explain the structure of dwarf galaxies [JO], we use a 
conservative bound [U] (see also [SJ [7] ) 

a/fh^ < 0.1 cm 2 /g (85) 

for characteristic velocities of 10 km/s. As we mention below, it would not be difficult to find parameters that violate 
this bound. 

In the usual treatment of dark matter, constraints such as this one, obtained from present-day observations, can be 
directly applied to bounds on physics at freeze-out or before. It is important to remember here that, in our model, 
the evolution of the chameleon field means that such a connection is far less direct, and such bounds typically do not 
apply in the early universe. 



D. Dark Decays 



The dark Higgs h and the dark gauge boson A are allowed to decay. As mentioned earlier, we assume m^(</) c ) > 
2Ma{4>c) so that h has a tree-level decay channel to A, as shown in Fig. [3} Its decay width is then 

r = P(<t>c) mg(^ e ) / 

Although the A particle is allowed to decay to r\ particles, which are substantially smaller in mass, this occurs through 
a 1-fermion-loop process, as shown in Fig. [4] The amplitude is also suppressed by two factors of m 2 from the 0(j] 2 ) 
term in the expansion of m^(<f>). The non-zero amplitude in the limit of m, <C Ma, rhip is 



4M|( 



™ 2 h(<f>c) 



12 m 



(86) 



14 







FIG. 4: 1-loop A decay. Only the first diagram is non-zero. 



4wr 



-k' 2 -e*{p) Afh^M A C [p 2 ,{p- k' 2 ) 2 ,kl,m^,m^,m^\ 



M = — 

+ (8to^ -I- M^)_B b 2 ,?Ti^,m^] - 8m 2 S [fc2 2 ,m^,m^] 

where B and Co are scalar Passarino-Veltman functions [42 45] , defined via 

1 



B [p 2 ,m 2 ,m 2 } = [d 
in 2 J 

Co[p 2 ,{p - Pi) 2 iPl^-m^m, to] = 



1 



d 4 Z 



{I 2 + m 2 )[(l+p) 2 + m 2 } 
1 



(I 2 + m 2 ) [{I + p) 2 + to 2 ] [(Z + pi) 2 + to 2 ] 
The Co integral is finite and, in the approximation m^, 3> M4 3> m n , reduces to 

C [p 2 , (p - k' 2 ) 2 ,k 2 , , to ^ , to ^ ] 



1 



4to 2 



(87) 



(89) 



(90) 



The Bo integral diverges, so we cut off the loop-momentum integral at some large scale. Using to 3 for this purpose, 
since we will find it numerically to be the largest mass-suppression scale in our theory, we have 



B [p ,fh^,m^\ « B [k' 2 ,to^,to^] ss 2 In 



TO 3 



Putting everything together, the decay width of A is then given by 



6M A \ to 2 J 



In I ^) 



1 2 



(91) 



(92) 



The A bosons must decay efficiently enough not to contribute significantly to the energy density budget today. Though 
the decaying exponential makes meeting this criterion more difficult, there is still sufficient parameter space for which 
the A energy density does not pose a problem. 



VI. NUMERICAL SOLUTIONS 



While we have described a number of ways to understand the evolution of the fields analytically, including, for example, 
the adiabatic approximation in which the chameleon tracks the minimum of its effective potential, ultimately, we are 
able to numerically solve the relevant equations of motion completely. To do so, of course, we must make sensible 
choices for our parameters to satisfy the various bounds and simplifying inequalities we have specified. 

We need to implement the correct relationship between the dark sector temperature T<j and that in the photon sector 
T, which in turn requires us to correctly enumerate the massless degrees of freedom at the relevant scales. At the 
unification scale, all the dark particles (ip, A, h, tp) are relativistic, so gfg(t u ) = 8.5. Around the epoch of dark 
matter freeze-out, only i\> is non-relativistic, so gfs(tf) = 5- Thus, at freeze-out, £/ = 1.19 for g*s(tf) — 106.75 or 
= 0.56 for g*s{tf) — 10.75. With these numbers, the bound on the number of effective neutrino species in Q is 
easily satisfied. 



15 



lO" 6 

10- 7 

lO" 8 

10- 9 

" lO" 10 

> 

a lo- 11 

lO" 12 
lO" 13 



10" 



10" 



10 



0.1 



0.001 



T [GeV] 



lO" 7 



10" 



lO" 11 



10" 23 
lO" 24 
lO" 25 

1Q -26 

lO" 27 £ 
lO" 28 & 
lO" 29 

1Q -30 



10 



-31 



FIG. 5: The evolution of the dark matter annihilation cross-section (av) as a function of T in GeV for the model 
with m^p = 125.4 GeV. The red, dotted line indicates the approximate dark matter freeze-out temperature. 
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FIG. 6: The evolution of the chameleon field as a function of T in GeV with = 125.4 GeV. The red, dotted 
line indicates the approximate dark matter freeze-out temperature. 



It turns out that our model is rather insensitive to Mj\((f)) and rh^cf)) at lowest order, and so we simply choose M^^q) 
to satisfy scattering cross section bounds. It is simplest to assume the attractor solution for <f> and then later verify 
that it is in fact adhered to. Interestingly, lower dark matter masses are more easily able to give large boosts to the 
annihilation cross section. Masses around 0.1 GeV < < 100 GeV can give boosts up to 4 to 6 orders of magnitude, 
whereas masses > 100 GeV give boosts that tend to be limited to 2 or 3 orders of magnitude. 

Finally, we must ensure that the evolution ends with the observed density of dark matter today. For this figure we 
use the bounds from the 7- year WMAP data j^Hj: assuming a ACDM cosmology 

Sl DM /i 2 = — = 0.1109 ± 0.0056 . (93) 

PcO 



Given these constraints, we numerically solve the Boltzmann equation subject to the parameter choices: = 
125.4 GeV; m T = m^; m 1 = 40 TeV; m 2 = 470 TeV; m 3 = 490 TeV; A = 20 GeV; A 2 = 0.7; A 3 = 9.8; g = 0.82; 
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FIG. 7: The evolution of the dark matter mass m,/,(</>) as a function of T in GeV with — 125 .4 GeV. The red, 
dotted line indicates the approximate dark matter freeze-out temperature. 
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FIG. 8: The evolution of the gauge coupling f(<f>) as a function of T in GeV with = 125.4 GeV. The red, dotted 
line indicates the approximate dark matter freeze-out temperature. 



and v = 8.2 GeV. This comprises an optimistic set of parameter choices, that satisfies all our bounds and provides 
a large change of order ~ 10 6 in cross-section over the history of the universe. The results for the evolution of the 
cross-section, the chameleon field, the dark matter mass, and the coupling / as a function of T — ttit/x are shown in 
Figs. [5j |6j [7J and[8j respectively. We also show the scattering cross-section in Fig. [9] With these particular choice of 
parameters, the scattering cross-section is too small to have interesting astrophysically observable consequences. 

However, there are other choices of parameters that will still give a boost to the annihilation cross-section while 
yielding a larger scattering cross-section. For example, here are another set of parameters: = 4 GeV; = to^; 
m 1 = 18 TeV; m 2 = 115 TeV; m 3 = 488 TeV; A = 83.5 GeV; A 2 = 8.6; A 3 = 7.5; g = 0.32; and v = 0.1 GeV. 
Again, all constraints are s atisfi ed, and the annihilation cross-section obtains a boost of ~ 10 between fre eze-out and 



today and is shown in Fig. 10 The scattering cross-section today is 0.37 cm 2 /g and is shown in Fig. 11 This value 
for the cross-section is the right magnitude to match observational bounds [7] • It is therefore possible to get a 
large scattering cross-section in the late universe while obtaining the right relic abundance in the early universe. This 
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FIG. 9: The evolution of the scattering cross-section a/fh^ as a function of T in GeV with m$ — 125.4 GeV. The 
red, dotted line indicates the approximate dark matter freeze-out temperature. 
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FIG. 10: The evolution of the dark matter annihilation cross-section (av) in a different model than the one plotted 
in previous figures. This one has m$ = 4 GeV and a larger scattering cross-section at late times. 

difference can ultimately be traced to the different ways in which f(<fi) and enter the annihilation and scattering 
cross-sections. 



VII. CONCLUSIONS 



In this paper we have investigated the possibility that the properties of dark matter depend crucially on the dynamics 
of a chameleon field - a scalar field whose cosmological evolution depends not only on its bare potential, but also on the 
local density of other matter (such as dark matter itself) in the universe. We have shown that such a coupling allows 
the annihilation cross-section (for example) of the dark matter particles to change by several orders of magnitude 
between freeze-out and today, while remaining consistent with all observational constraints. We have presented a 
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FIG. 11: The evolution of the scattering cross-section a/ih^ in the model with = 4 GeV. Note that while the 
final scattering cross-section is relatively large, it does not evolve significantly in the late universe. 



general formalism to describe how this might happen, and have provided a specific particle physics example, in 
which all relevant quantities can be calculated. While there are significant observational and theoretical constraints 
on models of this type, it is nevertheless possible for the cross-section to evolve in such a way that there may be 
interesting implications for the detection of dark matter, and for its dynamical effects on late-universe astrophysics. 

There are, of course, other possible complications to this idea that are beyond the scope of the current paper, but that 
provide interesting avenues for future study. One natural step is to couple our model directly to the Standard Model. 
One way to achieve this is to directly add the dark U(l) to the current SM gauge group [IB]. Another possibility is 
to couple to the Standard Model through U(l) kinetic mixing J37J EE]- This extension of our model should be able 
to easily accommodate the relevant particle physics constraints [49H52] , while allowing for decays of the dark gauge 
boson to Standard Model particles well before BBN. However, it is a more delicate issue to decide what a natural 
route would be to couple the visible and dark scalar sectors, particularly with regards to coupling the chameleon to 
normal matter. 

Finally, we did not attempt a careful analysis of the effect of late-universe inhomogeneities on the chameleon field or 
the dark matter properties on which it depends. In the specific models we considered, it seems as if such effects would 
be small, but a more careful examination is warranted. 
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Appendix: Feynman Rules 



The Feynman rules are shown in Fig. [12] All of these diagrams have higher order corrections that involve rj particles. 




\ 
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> = -«|A A v 
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/ 



FIG. 12: Feynman rules for h (dashed line), ip (solid line), (wavy line), and r\ (dotted line). We include the 
Yukawa interaction with ip and rj, which is relevant for the 1-loop A-decay amplitude in (87), but other 77-interaction 
vertices are not shown. All parameters labeled by a tilde are evaluated at <j) c . 
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